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Abstract 



o 



C ' In contrast to multigap superconductors (e.g. MgB2), the low-temperature properties 

of nodal superconductors are dominated by nodal excitations. Here we extend for a 
variety of nodal superocnductors the earlier work by Simon and Lee and Kiibert and 
Hirschfeld. The scaling relations seen in the thermodynamics and the thermal conductivity 
will provide an unequivocal test of nodal superconductivity. 

> 

CN \ 1 INTRODUCTION 

Although nodal superconductors have been with us since 1979 [1], the systematic study of the 
ly-j gap symmetry of these new superconductors began only around 1994 with the establishment of 
\ d-wave symmetry of high-Tc cuprate superconductors through the angle resolved photoemission 

■ spectrum (ARPES) [2] and Josephson interferometry [3,4]. Unfortunately, however, these 
powerful techniques have not been applied to other nodal superconductors like Sr2Ru04, heavy- 
fermion superconductors and organic superconductors. 

Since 2001 Izawa et al have succeeded in determining the gap functions |A(k)|'s in Sr2Ru04 
[5], CeCoIns [6], k-(ET)2Cu(NCS)2 [7], YNi2B2C [8], PrOs4Sbi2, [9,10], and UPd2Al3 [11,12] 
Q ' through measurements of the angle-dependent thermal conductivity in the vortex state. These 
O ■ experiments are only possible now since a) high-quality single crystals of these compounds are 
now available, b) low-temperature facilities which allow one to reach 1 - 0.1 K are available, 
and c)the necessary theoretical development following the seminal paper by Volovik [13]. 
. ^ Indeed Volovik's approach has been extended in a variety of directions, as reviewed in 

■ [14]. Also the angle dependent magentothermal conductivity and the scaling relations [15] 
■ ■ ' in the vortex state will provide a crucial test of nodal superconductivity. For example, the 

multigap superconductors do not exhibit the scaling relations we are going to discuss in general. 
Therefore, if any given superconductor exhibits a scaling relation discussed here, it is very likely 
that the material is a nodal superconductor. For example, the specific heat data of Sr2Ru04 
by Deguchi et al [15] obeys the scaling relation given in [16]. Therefore the simplest choice of 
gap function in Sr2Ru04 is the chiral f-wave superconductor as pointed out in [5]. 

The scaling relations in the vortex state in d-wave superconductors were first proposed by 
Simon and Lee [17]. Then within the semiclassical approximation, a la Volovik [13] Kiibert 
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and Hirschfcld (KH) [18] have succeeded in deriving the scahng function for the quasiparticle 
density of states. KH then calculated the thermal conductivity in the scaling region [19]. An 
error in [19] was pointed out and corrected in [20]. However, in [20] only the asymptotic 
behavior of the thermal conductivity (T <C< ]v • q] >, where v • q is the Doppler shift) has 
been worked out. 

In the following we shall derive the scaling relations for a class of quasi-2D superconduc- 
tors, where |A(k)| = A|/| and f = cos(2(/)), sin(2^) (d-wave superconductor), / = e^^^cosx 
(chiral f-wave superconductor as in Sr2Ru04), / = cos(2x) (g-wave superconductor as in 
UPd2Al3 [12].) These superconductors have the same quasiparticle density of states as in 
d-wave superconductors [21] 

N{E)/No = G{x) (1.1) 

where 

2t 

G{x) = —K{x) for X < 1 (1.2) 

TT 

= -K{x-'^) for X > 1. (1.3) 

TT 

where x = \E\/A and K(k) is the complete elliptic integral of the second kind. In particular 
for \E\ < 0.3A we have G{E/A) = Jfi. 

As discussed elsewhere [14], all the nodal superconductors so far discovered have G(E) 
~ \E\/A. Then one can establish a variety of scaling relations in the superclean limit [20], that 
is, for (FA)^/^ < T, i?; V • q < A where T, E, A and F are the temperature, the quasiparticle 
energy, the maximal value of the energy gap and the quasiparticle scattering rate respectively. 
Therefore the scaling relations provide another test for nodal superconductivity. 



2 QUASIPARTICLE DENSITY OF STATES 

Let us limit ourselves to a class of quasi-2D systems with f listed in the preceding section. As 
already noted we have G{E) ~ \E\/A for E <^ A. In the presence of a magnetic field we find 

G(£;,H) = <|£;-vq|>A-^ (2.1) 

where v • q is the Doppler shift and (...) means the average over the Fermi surface and the 
vortex lattice. When H || c in the class of quasi-2D systems, the average can be performed 
analytically and we find [18] 

G{E,U) = ^jr9{E/e) (2.2) 

where 

g{s) = J5(l + ^),s>l (2.3) 

= ^y/l- + ^(1 + 2s2)arcsin(s), s<l (2.4) 
4 4s 

where e = \v^/ eH and v is the Fermi velocity within the ab plane. The scaling function 
^G{E,'R) is shown in Fig. 1. As is readily seen G{E,'R) for H || c cannot discriminate 
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between different |A(k)|'s in the above class of nodal superconductors. Then the scaling 
function for the specific heat is given by [16] 

C7,(T,H)/C,(T,0) = F{T/e) (2.5) 

where 

.[l + 4(-r + T4(-)^ + -]>fore/T>l (2.8) 



97rC(3)T' 18' e ' 1800' e 



The scaling function and the experimental data for Sr2Ru04 [15] are shown in Fig. 2. As 
is seen readily the scaling function gives an excellent description of the experimental data. 
As noted in [15], this clearly shows the superconductivity in Sr2Ru04 is consistent with the 
chiral f-wave superconductor as discussed in [22]. On the other hand, as noted in [23], this is 
incompatible with p-wave superconductivity. 

Now when H is rotated within the a-b plane with an angle </) from the a axis, we can 
discriminate between / = cos(2(/)) and / = sin(20) (the case of vertical nodes). We obtain for 
/ = cos(2.^) 



G(^,H) = 




Prl-4 



JOURNAL DE PHYSIQUE IV 



10 - 



5 ■ 



a. 

y L 

c 

o 

O 1 L.. 



-Hx) 



J — I — I, I ,1 I r I 



4 B 8 2 

0.1 



-J 1 — I I I 1 1 > 



4 B S 2 4 B 8 

1 10 



Fig. 2. The scaling function F{T/e) — F{x) and Sr2Ru04 specific heat data from Ref. [15]. 




^(0.963 + 0.0205 cos(40) + i(1.132 - 0.081 cos(4(/))) x (^^^^ 




(2.11) 



(2.12) 



where e = ^vy/eH and w = y'vvc- For f=sin(20) we have the same forniulas as in Eqs. (12) 
and (13), except that cos(40) in Eq.(14) should be changed to — cos(40). Also the presence 
of the fourfold term in the specific heat has been studied by Revaz et al [24] . They found no 
fourfold term within an accuracy of 3%. This suggests strongly that the thermal conductivity 
provides a more sensitive test of the gap symmetry. 

For superconductors with horizontal nodes (e.g. f — sinx,cos(2x),cosx) the field configu- 
ration H [| b — c plane, with 9 the angle H makes from the c-axis, is more appropriate. Then 
we find [12] 



-E 



e±(^^,0,x) 



G{ 



E 



a(1+2^)'^-E<^ 



(2.13) 



(2.14) 



'l±]-sH24>)E{- 




± i sin(20) 



I ± i sin(2(/)) ^ 



(2.15) 
(2.16) 
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^(0.963 + 0.0205 cos(#) + ^(1.132- 0.081 cos(40))) x (-)^ + ...for ^ > (12.17) 
ttA 6 e E 

4e ^— 1 



-^^/x/2(l - sin^ 61 (sin^ (6*) + 16a^ cos^ 61 sin^ Xo)a;"^ + 
ttA Id 

i(-)2a;-i(l + 4sin2^(sin2(0) + 16a2cos2^sin2xo)a;-'))for | > 1 (2.18) 
3 e 16 hi 

where e = vVeH, a = Vc/v and x = 1 + cos^ ^ + 2a^ sin^ ^ sin^ xo and Xo = 0, ^ and ^ for 
/ = sinx,cos(2x) and cosx respectively. Therefore in the present configuration the angle 
dependent thermal conductivity can discriminate different A(k)'s with horizontal nodes. 

So far we have completely ignored the effect of impurity scattering. As already indicated 
the present analysis is valid in the superclean limit [14,20], i.e. for (FA)^/^ < T, i5, v • q < A 
where F is the quasiparticle scattering rate in the normal state. Then the superclean limit 
appears to require F/A < 0.01. 

3 THERMAL CONDUCTIVITY 

In the past few years the angle dependent magnetothermal conductivity (ADMTC) has proven 
itself the most powerful technique to probe the nodal structure of the gap function A(k). Also 
in many cases the nodal structure of A(k) is sufScient to deduce A(k) itself. We are concerned 
that much of the confusion and the controversy in the literature regarding the gap functions 
in Sr2Ru04, PrOs4Sbi2 and UPd2Al3 may be largely due to a misunderstanding of Volovik's 
approach. References [14, 20] contain a detailed description of this approach. Generalizing the 
standard expression of the thermal conductivity given in [25,26], the thermal conductivity of 
the class of nodal superconductors in the vortex state is given by [14] 

n r . .2/ M^,H) \_^2. 

where 



and 



2 V |(w-v-q)2- A2/^| 



r = Imv'(w - V • q)2 - A2/2 (3.3) 



Here < ... > denotes the averages over the Fermi surface and vortex lattice [20]. In the 
superclean limit a) is given by 

^ = o^ + zf/ ^ '"-"-q' \ (3.4) 

\V('^-vq)2-A2/2/ 

~ Lu + irG{uj,B.) (3.5) 
in the Born limit. And in the unitary limit we find 

u) = uj + irG-\LO,H) (3.6) 
where G(a;,H) has been defined in Eq.(5). 



Prl-6 



JOURNAL DE PHYSIQUE IV 



First we limit ourselves to the quasi-2D systems in a magnetic field H || c. The in the Born 
limit we obtain 



f°° , 2 , 2/ w T^^nT 1 
K = - — / au CO sech (— - = — — — = -k„ (3.7) 
SmT^T Jq ^2T 12mr 2 ^ ' 



where k„ = g^^p is the thermal conductivity in the normal state. In particular Eq.(28) gives 
the scaling function 

The last result agrees with the corresponding result given in Rcf. [19] despite the use of a rather 
unphysical spatial average in this work. In the unitary limit, on the other hand, we obtain 



8m(TA)2r 
n 



POO 

/ cL;a;^ < Iw- v-q| >^ sech^(a;/2T) 



(3.9) 



(-^)^ / dwu^G^{u/e)sech^{u/2T) (3.10) 



8mT2r'4A' 

where G{w/e) has already been defined in Eqs. 5 and 6. This has asymptotics 



60mrA2 \^ 7'ttT' 28^ttT' 

TT^nT , 7re , 9 777^ T „ , /„,„s 

(tx) 1 + ^(t) +•■• ' fore»T (3.12) 



12TOr'4A' V 15 ' e 
where e = 2iv^2 . Then the scaling function is given by 



f 1 + + . . .V fore » T (3.15) 



112'A' V 15 

These scaling functions are shown in Fig. 3. In this figure we also include the scaling function 
when the inversion symmetry is broken in the impurity scattering [27]. Fi(T/e) describes the 
thermal conductivity data of the non-centrosymmetric triplet superconductor CePtsSi by Izawa 
et al [28]. The scaling function F„(T/e) is very different from the one given in Ref. [19] but 
describes consistently the scaling behaviors of the thermal conductivity of UPts as reported 
by Suderow et al [29]. 

4 ANGLE DEPENDENT THERMAL CONDUCTIVITY TENSOR 

Let us consider dj.2_j^2-wave supcrconducivity as in the high-Tc cuprates, CeCoIns [6] and 
k;-(ET)2(NCS)2 [7] in a magnetic field within the a-b plane. Then the thermal conductivity 
tensors within the ab-plane are given by 

= K,, = g^^°°da;a;2(^^^^G(a;/e(0,x))^'sech2(a;/2r) (4.1) 
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where e(0, x) = sin(2^) — | cos(2x))^^^ and (. . .) means the average over ± and over 

X- Here wc assumed the unitary impurity scattering and the superclean hmit in the present 
derivation. This gives the following asymptotics: 

7nT ,7rr.n 5, e ,9 15 , e \ „ /. „s 



AnT e 7 ttT 

= ^^(—)2[0.927 + 0.039 cos(4(/)) + —( — )2[1.090 - 0.055 cos(4(/.)]for e > 1(4.3) 

First of all, the present result is consistent with that in Rcf. [20] for T < e. On the other hand, 
for T > e there is no fourfold term. In other words the present theory in the superclean limit 
cannot describethe fourfold symmetry in Kxx observed in YBCO for T > UK [30-32]. We 
have proposed the sign inversion of the fourfold term for T > e in the clean limit in [33] . 
Also the Hall conductivity in the present geometry is given by [20] 

/•OO 

/ (L;a;^(sin(2^')|a;-v-q|)(|a;-v-q|)sech2(a;/2T) (4.4) 
Jo 



8m(TA)2r 
Further, we find 

(sin(20)|a;-vq|) = - sin(20)— -, for e «; [wj (4.5) 

2\uj\ 

~ -sin(2(/))^(0.535- (-)20.14192)for e» |w| (4.6) 
Inserting these into Eq.(40) we find 

Kxy = -sin(2</.)24^(^)^(l+2(;^)',fore«T (4.7) 

AnT e 7 ttT 

= -(—)2sin(20)(O.5152 + 0.011 cos(4<?i)) ( — )2 (0.213 + .0607 cos(4(?i)) fore (*.ffl 

3mr A 30 e 
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Therefore in the superclean limit K^y ^ — sin(2(^)iJ independent of e/T. Also as ^ increases 
the coefficient of -sm{2(j))H decreases almost 40%. The present result appears to be consistent 
with the Hall conductivity data of YBCO reported by Ocaha and Esquinazi [34] for e/A < 1. 
Also in the superclean limit the sign of the Hall conductivity is the same for all T as long as 
T < A(T). 

5 CONCLUDING REMARKS 

We have shown a) the thermal conductivity in nodal superconductors for T < O.STc is domi- 
nated by the quasiparticles or nodal excitations, b) the quasiparticles in the vortex state are 
accurately described in terms of the semiclassical approximation. Thus the angle dependent 
magneto-thermal conductivity provides a powerful tool to determine the nodal structure of the 
gap function as demonstrated by a series of experiments by Izawa et al [5-12]. Also in most 
cases the nodal structure of the gap function is adequate to deduce |A(k)| itself. In addition 
we have shown that all these model superconductors exhibit a variety of scaling relations. We 
have proposed scaling relations for PrOs4Sbi2 [35]. Furthermore, from the unusual scaling 
relation seen in the thermal conductivity in CcPtaSi we can deduce anomalous impurity scat- 
tering in this system lacking crystalline inversion symmetry. Indeed the scaling relations in 
nodal superconductors provide a unique way to characterize this new class of superconductors. 
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